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1 Introduction 



The quasi-geostrophic(QG) equation models large scale geophysical flows. It is derived as 
an approximation of the rotating Navier-Stokes equations by an asymptotic expansion in a 
small Rossby number. The barotropic QG equation is written in terms of stream function 
^(x,y,t) (Q, 16, p. 234): 

+ J(V, AV) + pfa = ^A 2 V - rA^ + f(x, y, t), (1.1) 

where (3 > the meridional gradient of the Coriolis parameter, v > the viscous dissi- 
pation constant, r > the Ekman dissipation constant, f(x,y,t) the wind forcing, and 
J (/) 9) = fx9 y ~ f y 9x is the Jacobian operator. 



Equation (IT) can be rewritten in terms of the relative vorticity lu(x, y, t) = Aip(x, y, t) 

as 

uj t + J (t/j, lo) + f3ip x = vI\lo -ru + f(x, y, t). (1.2) 

For an arbitrary bounded planar domain D with sufficiently regular (such as, piecewise 
smooth) boundary dD this equation can be supplemented by homogenous Dirichlet bound- 
ary conditions for both ip and u, namely, the no-penetration and slip boundary conditions 
proposed by Pedlosky ||, p. 34: 

ip(x,y,t)=0, uj(x, y, t) = 0, on dD, (1-3) 

together with an appropriate initial condition, 

u(x,y,0) =u (x,y), inD. (1.4) 

The global well-posedness (i.e., existence and uniqueness of smooth solution) of the 
dissipative model (1.2)-( |l~4| ) can be obtained similarly as in, for example, ||, ||], ||; see 



also ||. Steady wind forcing has been used in numerical simulations (0]). Brannan et al 
[s] considered the effect of quasi-geostrophic dynamics under random forcing. Duan et al 



U and [JTOj] obtained the existence of time periodic, time almost periodic quasi-geostrophic 
response of time periodic and time almost periodic wind forcing respectively. 

In this paper, we assume that the right-hand side or the forcing term of the QG flow 

model ( |1.1[ ) is rapidly oscillating, i.e., it has the form f(x,y,t) = f(x,y,rjt) = f(r]t), 
with parameter > 1. We also assume that / has a time average. With such forcing, it 
is desirable to understand the fluid dynamics in some averaged sense, and compare the 
averaged flows with the original unaveraged flows. 

The main result of this paper is the averaging principle for quasi-geostrophic motions 
with rapidly oscillating forcing, both on finite but large time intervals and on the entire 
time axis. This includes comparison estimate, stability estimate, and convergence result 
(as rj — > oo) between quasi-geostrophic motions and its averaged motions. We also inves- 
tigate the existence of almost periodic quasi-geostrophic motions under almost periodic 
forcing, and the convergence of the attractor of the non-autonomous equation (|1.1| ) to the 
attractor of the averaged autonomous equation as r\ — > oo. 
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In §2, we study the averaging principle for the QG flow model on finite but large time 
intervals and in §3, we extend the averaging principle to the entire time axis. In the rest of 
this section, we briefly review some background and provide some preliminaries for later 
use. 

Starting from the fundamental work of Bogolyubov [11] the averaging theory for ODE 
has been developed and generalized in a large number of works (see [ 12 1 — 114] and the 
references therein). Bogolyubov's main theorems have been generalized in Jl5| to the 
case of differential equations with bounded operator-valued coefficients. Some problems of 
averaging of differential equations with unbounded operator-valued coefficients have been 
considered in |16||-[p^] in the framework of abstract parabolic equations. In | |19[ ] , Ilyin 
considered the averaging principle for an equation of the form 

d t u = N{u) + f{ V t), (1.5) 

where / is a given right-hand side and r] 1 is a large dimensionless parameter, and / 
satisfies 

lim \ f f(s)ds = f . (1.6) 

t^co t Jo 

Standard abbreviations L 2 = L 2 (D), Hq = Hq(D), k = 1, 2..., are used for the common 
Sobolev spaces, with < •, • > and || • || denoting the usual scalar product and norm, 
respectively, in L 2 . We need the following properties and estimates (see Q) of the Jacobian 
operator J : Hq x Hq — > L 1 : 

/ J(f,g)hdxdy = - J(f,h)gdxdy, J(f,g)gdxdy = 0, (1.7) 

JD JD JD 

| / J(f,gdxdy\ < ||V/||||Vg||, (1.8) 

JD 

for all f,g,h£ Hq, and 



| J D J(Vf ig )Ahdxdy\ < WJ-I||A/||||A<7||||A/»||, (1.9) 

for all f,g,h € Hq. We also recall the Poincare inequality pH] 

\\ g f < H f \Vg\dxdy = ^||V<?||, (1.10) 
vr Jd vr 

for g G Hq, and the Young's inequality p0[| 

AB < -A 2 + — B 2 , (1.11) 
~ 2 2e v ' 

where A, B are nonnegative numbers and e > 0. 

We can further rewrite the QG flow model (|1.2|). From 



A^ = uj, (x,y)£D, <p\ dD = 0, (1.12) 
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we get tp = A 1 lo. Thus fll-2| ) can be rewritten as 

oj t + J(A _1 u;, w) + (3d x l\~ 1 ijj = i/Aw - ru + /(x, y, i). (1.13) 

Let 

-A = vA-rI-0d x A~ 1 . 

Then by a result in |fl6|1 , we know that A is a sectorial operator, and hence e - " 4 * generates 
an analytic semigroup in I? . 

We will give a sufficient condition to ensure the smallest eigenvalue of A to be positive. 
Consider the eigenvalue equation Au = Xu. We have the following energy estimate 



\\\u\ 1 2 = u\\ Vull 2 + Hlull 2 



/ A 1 ud x udxdy 
Jd 



i d 

> z;|M| 2 + r||-u|| 2 - /3\\ A _1 u|| Vu|| 

. ||V7 ,12 , || || 2 P\D\ „ I, 

> v vu + r \\u\\ a Va 

7T 

\ IIV7 1 1 2 , II ||2 , I, ,|2 . Mv7 ||2\ 

> i/||Vu|| + r it ( a ilMI + a 2||Vit|| ) 

7T 

> („ - PM ai )\\Vuf + (r - ^a 2 )|M| 2 , (1.14) 

7T 7T 

where we used the Poincare inequality ( 1.1 Q| ) after the second inequality sign above, and 
where arbitrary constants ai,a 2 satisfy aia 2 = \. Therefore, when 

B 2 \D\ 2 

Aur > ^-V". ( L15 ) 



7T 2 



and if we take a 2 = then we have 



P 2 \D\ 2 , n „_ ll2 ^ 2 \D\\\D\ n 



A|HI - ( " " ^ )I|Vm|1 - ~v M ■ (L16) 

So, when z/, r, /3 and \D\ satisfy the condition ( [L.15 ), the smallest eigenvalue of A is positive. 
In this case, the QG flow model 

u>t + Auj + J( A _1 w, uj) = f(x, y, t) 

is a dissipative dynamical system. 

We note that the condition ( 1.1EQ is sharper than the corresponding condition in Duan 
and Kloeden |L0| ] . We define the fractional power of A as follows |lq] : 

A a = (A~ a )~\ where A~ a = — — / t^e^dt. 

r(a) Jo 

The corresponding domains D{A a ) are Banach spaces with norm given by 

\\x\\ a ■= W^Hd^A") = W-^xW. 
For the rest of this section, we recall some definitions and useful results for later use. 
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Theorem 1.1 Jig / The following estimates are valid: 

\\e~ At \\ L 2^ L 2 < Ke~ at , t>0, (1.17) 

\\A a e- At \\ L ^ L2 < ^e~ at , t>0, (1.18) 
where K, K a are positive constants. 



Theorem 1.2 J?6| / Given two sectorial operators A and B in I? , let D(A) = D(B), 
Rea(A) > 0,Rea(B) > 0, and for some a € [0,1). Let the operator (A - B)A~ a be 
bounded in L 2 . Then for every 7 € [0, 1), D{A y ) = D^B 1 ), the two norms being equivalent. 



Definition 1.3 [VY\ A continuous function f : IR — > X is called almost periodicfa.p.), 
if for every e > there exists a number I = /(e) > such that each interval (T,T + /) 
contains a point r = r e (called an almost period) satisfying the inequality 

Su PtelR \\f(t + r)-f(t)\\<e. 

If / depends on other arguments, then the above inequality holds uniformly with respect 
to norms. 

It follows from the theory of a.p. functions that there exists a countable set of number 
A Q for which 



-r 



hm -L f fity^dt + 0. 

T— >oo Zl J—T 



The numbers {X a } are called the Fourier exponents of /|17|. 



Definition 1.4 J7j/ A countable set of numbers {ui a } is called the frequency basis (denoted 
by Mf) of an a. p. function f if every X a can be uniquely respected by a linear combination 
of the numbers uj a with integer coefficients. 



Definition 1.5 [VY\ For a given a.p. function f , the sequence {t m } is called f— current 
if 

^PtelRll/^ + *m) - /(*)|| < £m -»• 0, as m ^ 00. 



Theorem 1.6 [l r /j Given two almost periodic functions f and g, suppose that every 
f— current sequence is also a g— current sequence. Then the frequency basis of g is con- 
tained in that of f: A4 g C Mf- 

We now turn to the averaging principle for the QG flow model. 
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2 Averaging Principle on Finite Time Intervals 



In this section, we consider the averaging principle for the QG flow model on finite (but 
large) time intervals. In the next section, we extend the result to the entire time axis. We 
assume that the right-hand side or the forcing term of the QG flow model ( [LIT ) or ( |1.2[ ) 
is rapidly oscillating, i.e., it has the form f(x,y,t) = f(x,y,rjt) = f(r/t), with parameter 
rj ^> 1. Let ?] > 1 be a large dimensionless parameter. Setting 

T = rjt, e = rf l , 

we obtain the equation in the so-called standard form 

lo t + eAtv + eJ(A _1 tj, u) = ef(x, y, t). (2.1) 

We assume that / has a time average in D{A 1 ); the value of 7 will be specified later on. 
More precisely, let /(r),/o S A' 7 and suppose that 

\\A\^f t +T f{T)dr-f )\\ <min(M 7 ,cT 7 (T)), (2.2) 

where M 7 > 0, cx 7 (T) ^ 0, as T — > 00. 
We consider the averaged equation 

u> T + eAu) + eJ(A~ 1 u), Co) = efo(x, y). (2-3) 



By the method and result of [21], |22[ and fl23f , we obtain the semigroup St corresponding 
to equation (|2.3|) possesses absorbing sets in the space H = L 2 ,V = D(A2) = Hq and 
D(A), || • || and || • ||d(^) denote the norm in V and D(A). These sets are certain balls 
B(R ) in these spaces, where Rq is large enough. This means that for every bounded set 
B 

S t BcB{R ), fort >t {B,R ). 

In addition, the semigroup is uniformly bounded in these spaces, that is, given any ball, 
in particular, the ball B(R$), there exists a ball B(R) such that 

S t B(R ) c B(R), fort > 0. 

By increasing R we may assume that 

S t B(R ) C B{R - p), for t > 0, p > 0, 

where p is a positive constant. We first consider the averaging principle in the space V. 
Given a point ujq in By(Ro), we compare the trajectories (solutions) cj(t), Cv(t) of system 
(2.1) and (2.3) starting from this initial point. Consider their difference on the interval 
t E [0, y]> T being arbitrary but fixed. We suppose for the moment that uj(t) £ By(R). 
Then the difference z(t) = uj(t) — Cj(t) satisfies the equation 

d t z + eA + e[J(A-^, a;) - J(A _1 a), Co)] = e(/(r) - / (r)). (2.4) 

We first have some estimates on the nonlinear terms. 
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Lemma 2.1 The nonlinear operator J(u,v) is a bounded Lipschitz map in the following 
sense: 

|| J(ui,vi) - J(u 2 ,V 2 )\\ < 
Ci(|Mi + \\u 2 \\x + \\vih_ + IMiXHui - + ||«2 - v 2 \\i_), (2.5) 

2 2 2 2 2 2 2 

\\J{ui,v t ) - J(u 2 ,v 2 )\\i < 

2 

Co(||«i||d(^) + IKHd(.A) + \\vi\\d(A) + INId(.a))(IN ~ v i\\d(A) + ll«2 - v 2 \\ D{A) ), (2.6) 
where Ci and Cq are some positive constants. 

2 

PROOF. Since 

J(ui,u 2 ) - J(vi,v 2 ) = 

{U\X - Vlx)U 2 y + (Uly - Vly)Vl x ~ (ll 2x ~ V 2x ) U\y ~ (u 2 y ~ V 2 y ) V 2x , (2.7) 

(2.5) and (2.6) are obtained by direct estimtes. Here the equivalence of norms || • \\jj2 and 

II • \\d(A) is used - ■ 
Now we get back to equation ( |2.4| ). Inverting the linear operator we come to an 
equivalent integral equation 



z(t) = e e- eA ^- s) [J(A- 1 u,Lo)-J(A- 1 LO,u;)]ds 
Jo 

+ e [ T e-^- s Hf( S )-f )ds. (2.8) 



Using (1.18) and (2.5), the D(^.2)-norm of the first term in the right hand side satisfies 
the inequality 

||e f A^e- eA{r - s \j{A~ l uj,uj) - J(A~ 1 uj,u})]ds\\ 



< e f T Kie-^(T - s)-h- ea ^ T - s hR\\z(s)\\ids 
Jo 2 2 

= 2RKie? [ T (t - s)-h- ea ^- s h\z(s)\\ids. (2.9) 
2 Jo 2 

Let us estimate the second term in the right hand side of (2.8). Integrating by parts we 
have 

"\-^(T-S ){f{s) _ fo)dsh 



ee -eA(r-s) {m _ fo)dtl r + e 2 Ae ^A(r-s) ( /(s ) _ /o ) ds || , 
J s J J s 

<\\eA^e-' M A^ j\f{t)-h)dt\\ 
Jo 

+ ||e 2 f T A^e-^-'^A 1 f T (f(s) - f )ds\\. (2.10) 

JO Js 
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Using (1.18) and (2.2), we further have 

WeA^e-^A* [ T (f(t)-f )dt\\ <eKi e- £aT (er)^\\- f A\f {t) - f )dt\\ 

JO 2 ' T Jo 

= (er)^Ki e~ £a T \\ - f A^f(t) - f )dt\\ 

27 T Jo 

< (er)H7Ki_^ m i n (Af 7 , C j 7 (T))e- ear =: L(r). (2.11) 

For any S > 0, let ts be so large that for r > Tg, <r 7 < 5. Let eo be so small that for e < eo 
then inequality — > Tg is valid. Then 



(eT)*+*Ki M 7 , if r < r 5 . 



L(r) <G 7l (T,e) = e- ea 

Let 7 > — \. Since does not depend on e, we let 8 — > and then e — > 0. We obtain 

llee - ^ [ {fit)- fo)dt\U < G 7 i (T, e) ^ when e ^ 0. (2.12) 
jo 2 

||e 2 f A^'e-^-^A 1 f T (f(s) -fo)ds\\ < Kzf^ 1 f mm{M 11 o-Ju))u'^du 

Jo Js 27 JO 

r — 

2 7 Jo 27 Jo 

= i^f _ 7 (7 + ^) _1 ((^)^ +7 + /^ l+7 ) = : GAT, e), (2.13) 

where for any [i > we have chosen so large that o~ 7 (t) < // when r > r„. Letting 
— > and then e — > we obtain 

G 72 (T,e)^0, e^O. 
Thus, by (2.8)-(2.13) we obtain the following inequality: 

||z(t)||i <Ke* f (r - 8)-*\\z(s)\\ids + GJT,e), (2.14) 
2 Jo 2 

where K = 2RKi and G 7 = G 7 i + G 72 -> 0, e -> 0. 
We need the following fact. 



Lemma 2.2 J7|/ Let 7 € (0, 1] and /or t € [0, T] 

<a + b [ (t- s)" / ~ 1 u(s)ds. 
Jo 

Then 

u{t) <oS 7 ((6r(7))7t), 
where the function E^(z) is monotone increasing and E^(z) ~ 7~ 1 e 2: as 2 — > 00. 
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Applying this lemma to the inequality (2.14) on r € [0, -j], we obtain 

\\z(t)\\i < G 7 (T,e)Ei(eTTrK 2 ) < G~ f (T,e)Ei(TTrK 2 ) =r&(e). (2.15) 

Using (1.18) and (2.6), we can do the same in D(A) assuming 7 > in (2.2) and obtain 

\\z(t)\\ D(A) < F 1 {T,e)E,(eT7rK 2 ) = 7#(e), F 7 (T, e) -> as e -> 0. (2.16) 

We thus have proved the proximity of solutions of (2.1) and (2.3) in V and -D(-4), assuming 
that the trajectory u(t) with initial condition u(0) G -By,D(.A)(-Ro) stays in the ball B(R) 
on the interval [0, — ]. 

Let e be so small that the right-hand side of (2.15) and (2.16) are less than ^, where p is 
defined earlier in this section when we discuss absorbibg sets. Suppose that the trajectory 
uj(t) leaves the ball B(R) during the interval [0, j] and let r* be the first momemt where 
||w(r*) = R. However, on the interval r G [0, t*\ both trajectories stay in the ball B(R) 
and what we have proved so far shows that the inequality ||w(r) — £>(r)|| < f is valid. In 
particular, it is valid for r = r*. This together with the inequality ||a>(r*)|| < R — p, which 
holds by the hypothesis of the following theorem and the property of the semigroup S(t), 
gives the contradiction 

Hr*)||<HO-w(0|| + ||w(0||<12-|. 

Therefore we have the following main result in this section. 

Theorem 2.3 (Averaging principle on finite time intevals) Let the right-hand side of 
equation (2.1) has an average in the sense of (2.2). Let T > be arbitrary and fixed. 

If j > — i and oj(0) = u)(0) G Bv(Rq), that is, the initial values coincide and belong 
to the absorbing ball, then for t € [0, j], 

\\lo(t) - Q(t)\\i < Vt(c) -» as e -» 0. 
7/7 > 0, and w(0) = cj(0) G B D ^(Ro), then for r G [0, j], 

||w(r) - w(r)|| D( ^) < r7r(e) -> as e -> 0, 
where r/-(e) and rjj,(e) are defined in (2.15), (2.16), respectively. 

This theorem gives comparison estimate and convergence result (as 77 — > 0) between 
the QG flows and averaged QG flows, on finite but large time intevals. 

3 Averaging Principle on the Entire Time Axis 

Now we turn to averaging principle for the QG flows on the entire time axis. Consider 

u> T + eAu> + eJ(A~ 1 cj, lo) = ef(x, y, t). (3-1) 
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All the hypotheses concerning the data of the problem are the same as those in §2; in 
particular, the average fo(x,y) exists in the sense of (2.2). 
We first consider stationary averaged QG flows: 

Au + J(A- 1 u ,u ) = f . (3.2) 

Note that this is not the stationary QG flow model, but the time-independent version of 
the averaged QG flow model. 

Under the condition of (1.15) and using Leray-Schauder fixed point theorem and el- 
liptic regularity ([pO|l), we know the stationary averaged QG flow model (|3.2|) has a unique 
stationary solution. Note that this unique stationary solution is denoted as Uo(x,y) and 
it is not to be confused with an initial datum for the time-dependent QG flow model. 

We change the dependent variable (from u> to z) in equation (3.1) via: 

U> = Uq + Z. 

Then by (3.2), we find that z satisfies the equation 

d t z = e(-Az - J(A- V uj) + J(A- V), wo) + fir) - f ). 

Since 

J(A~ x u),u) - J(A' 1 uq,u q ) = J(A _1 o;,z) + J{A~ x z,lj q ), 

we have 

8 t z = e(-Az - J(A~V z) - J^A-^uq) + f(r) - /„). 
Changing the dependent variable once again (from z to h) via, 

z = h — ev(r, e), 

we obtain 

d t z = d t h - ed t v = e(-Ah + eAv - J(A~ V h) - J(A" 1 /i, uj ) 

+ (eJ(A-^, v) + eJ(A- 1 v, uq) + f - fo). 
We chose the auxiliary function v(r,e) to satisfy the equation 

8 T v = -eAv + f -f. (3.3) 

Then we obtain the following equation for the new dependent variable h: 

3 T h = -e(Ah - J(A- 1 w , h) - J{A- l h, w )) + 

e(-J(A" 1 /i,/i) + J(ev,h) + eJ(A" 1 (w - ev),v) + eJ{A~ 1 h,v) + eJ(A" 1 w, o; ))- (3.4) 

For the rest of the section, we study the equation (|3.4|) for the new dependent variable 
h. Now, we first consider equation (3.3) for the auxiliary function v(r, e). 
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Lemma 3.1 Assume that the function f has an average in the sense of (2.2). If a— 7 < 1, 
then equation (3.3) has a unique solution v(r,e) bounded in D(A a ) uniformly in r € R. 
Moreover, 

||eu(r, e)|| a — > 0, as e — > 0. (3-5) 

If f is almost periodic with values in D{A~'), then v is almost periodic in D(A a ) with 
frequency basis contained in that of f . 

Proof. The desired solution is given by the formula 

v(r,e)= f T e -^- s \f -f))ds. (3.6) 



The uniqueness will be proved in a more general context in Lemma 3.2. below. Now we 
prove (3.5). Integrating by parts and using (2.2) and (1.18), we obtain 



kr(~Oi:„ IK / e' eA{r ' s) {fo-f))ds\ 

' —00 

fOO 



e' eAs (f -f(r-s))ds\\ a 

poo ps 

= ||e 2 / Ae~ eAs (f -f(T-t))dtds\\ a 
Jo Jo 

poo ps 

= ||e 2 / A^^'e-^sis^A 1 / (f - f(r - t))dt)ds\\ 
Jo Jo 

POD 

< e 1+a ^K l+a ^ / s^- a e- eas mm(M y ,a{s))ds 

Jo 

< K 1+Q _ 7 (e 1+a ^M 7 (l + a - jyh'f^ 1 + 5a 1+a ~"'T(l + a - 7)), 

where 5 is small and so = so(5) is so large that a(s)5 when s > sq. Letting 5 — > and 
then e — > 0, we obtain (3.5). 

Finally, let us prove the last statement of Lemma 3.1. By Theorem 1.6, it is sufficient 
to show that every /—recurrent sequence {r m } is also v— recurrent. By (3.6), 

POO 

V(T + Tm ) ~ V(T) = / e~ eAs (f(T -S)~ f(T + T m - S ))ds. 

Jo 

Therefore 

su PreR\\ v (. T + T m) ~ v( T )\\at < sup T6fl ||/(r + r m ) - /(r)|| 7 i^ a _ 7 r(7 - a + l)(ea) Q ~ 7 . 

Thus {r m } is indeed v— recurrent. The proof of Lemma 3.1 is complete. ■ 
Now we go back to study the equation (3^) for the new dependent variable h. We 
consider the operator 

Ch = Ah- J(A _1 /i, lo ) - JiA^cuo, h). (3.7) 
The operator C + XI has a compact inverse in H = I? for A > Ao = 2 a^|Z3| H^oll 2 ' wnere 

A -y 2 \D\ 2 

Ai = v — U— j ■ In fact, we consider the equation 

Ah - J{A~ l h, cj ) - ^(A _1 u;o ; h) + Xh = f. 
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Multiplying by h in I? and observing that (J(A o>q, h),h) = 0, we obtain 

(Ah, h) - (J(A _1 /i, wo), /i) + A||/i|| 2 = (/, h). 
Using ff 2 we find that 



^l|v/ i || 2 + ^ 

2 2 2Ai 

By (1.16) and (3.2) we have 

7T 

For A > Aq, ((£ + X)h, h) is coercive and by the Lax-Milgram lemma the equation 



-(J(A-%u ),h) < ||V/i||||Vwo||iM < -±\\Vhf + — \\h\\ 2 \\Vu \\ 2 . 



(Ah,h) > Xi\\Vh\f, HVwoll 2 < -jj— ||/o 



((C + X)h,h) = f 

has a unique solution h G D(.4.2). Note that the embedding D(A?) — > L 2 is compact. 
Hence C has a compact resolvent. Using (1.9), we can estimate J(A _1 /i, cjo) as follows. 

|J(A -1 /i,u;o) + J(A _1 a;o,MI < INI ll Vw o|| + ||w ||||V/i||. 

The operator A is sectorial, and due to the above estimates , we see that the operator 
C is also sectorial ( ]l6|). Moreover, if ReX > Xq, then the operator C\ = C + XI is 
invertible; hence 

Rea(C Xo ) >0. 

We also noe that 

(£ Ao - A)A~ a = X a A~ a - JiA^A- 1 -,^) - J(A- 1 w ,^" a -)- 

If a > 0, the operator £a is bounded from 1? — ► L 2 (using (1.9)). The fact that -D(^4) = 
D(C\ ) follows from the theorem on the regularity of the 2-order elliptic problem[(pO|], 
It now follows from Theorem 1.2 that D(A a ) = D(£^),a G [0, 1], that is 

ci a \\C a Xo \\ < \\A a \\ < c 2a \\Ct \\, V/i G £>M Q ). (3.8) 

We have thus shown that £ has a discrete spectrum. Let the stationary solution lvq be 
such that Re(£) ^ (which depends on the choice of /o). In other words, we suppose 
that 

a{£) = <t+(£)U<t-(£), o + (C)C\cr-{C) = and Rea + (£) < -a,Rea-{£) > a, a > 0. 

We observe that <x+(£) is a finite set of eigenvalues. The + sign indicates the unstable 
modes and — sign indicates the stable modes. 

Let 7_|_ be a coutour in the left half-plane enclosing a + (£). We set 



p ^ = iiIj XI - A >" dX - p - =1 - p * 
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The operator C can be decomposed as 

£ = £+ + £_, £+ = P+C, £- = P-C, 

where P + C C CP+, P^C C £P_. 

Setting if + = P+H, dimH+ = iV < oo and iT_ = P_H, we see that 

£+#+ C C 

The operator £ + G C(H + ) (bounded linear operator space on H+) and — £_ generates 
an analytic semigroup in iP_. Note that P+ and P_ commute with £ and the semigroup 
T(t) (which is generated by — C), in the sense that P+C C CP + ,P~C C CP-, P + T(t) = 
T(t)P+ and Pr(t) = T(t)P- for t > 0(||]). 

We consider the following equation for i G IR: 

d t h + Ch = fit). (3.9) 



Lemma 3.2 J7|/ Let /(t) G L^IR; D(/7)), 7 > 0. //a - 7 < 1, tfien £/ie equation (3.9) 
has a unique solution h bounded in D(C a ): 

\\ h \\c b QR;D(£<*)) < K ( a n)\\f\\ Lx (\R;D(^))- ( 3 ' 10 ) 



We continue to study equation (3.4) . Let F(h,e,r) = — J(A 1 h,h) + J(e v,/i) + 
eJ(A~ 1 (wo — e + eJ(A~ 1 h,v) + eJ(A _1 f , uq), then for /ij G D(C^) and assuming 

that < p(i = 1,2), we have 



||F(fc lj e,T)-F(ft 2j e,T)|| < iVi (e,p)||/»i - h 2 \\i, (3.11) 

2 2 

\\F(0,e,r)\\<m(e), (3.12) 

where iVi (e, p), Ni (e) — ► as e, p — > 0. 

2 2 

The proof of (3.11) and (3.12) can be obtained using Lemma 2.1 and Lemma 3.1. 
Moreover, if hi G D(C)(i = 1,2), using Lemma 2.1 and Lemma 3.1, we can prove 
F : D(A) — ► D(A^) is a bounded Lipschitz map (using Theorem 1.2). 
Converting equation (3A) to the original time variable t = er, we obtain 

dth + Ch = Q(h,e,t), (3.13) 

where Q(h,e,t) = F{h,e, ~). Obviously, Q(h,e,t) satisfies (3.11) and (3.12). 

Lemma 3.3 Assume that Q satisfies (3.11) and (3.12). Then if e < eo for some €q > 
small enough, (3.13) has a unique bounded solution h* with the following properties: 
1. \\h*\\ , 1 < 6(e) ^0, ase^O. 

11 n C b QR;D(A1) - W ' 
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2. There exists an initial manifold M~ in D(A 2 ), codimA^- = N, such that if 
ho E M- and \\ho\\i < p for p small enough, then the solution h(t) of equation (3.13) 

2 

with h(0) = ho satisfies the estimate 

\\h(t) - h*(t)\\i < Ke- aW \\h - h*(0)\\i, a > 0, as t -» oo. (3.14) 

2 2 

In particular, if N = 0, then the solution h* is asymptotically stable. There exists an 
initial manifold M. + in D(A^), dimj\4 + = N, such that (3.14) holds as t — > — oo. 

3. If the function Q(h,e,-) : IR — > D(A?) is almost periodic, then the solution h* is 
almost periodic with frequency basis contained in that of Q. 

This lemma is a special case of |19| |, Its proof is based on the contraction map principle, 
stability argument, Lemma 3.2 and Theorem 1.6. 

We now have the following main result in this section. 

Theorem 3.4 (Averaging principle on the entire time axis) Assume that the forcing f on 
the right-hand side of the quasi- geostrophic flow model (2.1) has an average in the sense 
of (2.2) uniformly in t € IR. Assume also that the spectrum of the linear operator C in 
( 3. 7) does not intersect with the imaginary axis. Then for e < cq small enough: 

1. In a small neighbourhood of the stationary averaged quasi- geostrophic flow ujq, the 
full quasi- geostrophic flow model (2.1) has a unique solution lo*(t), which is bounded on 
the entire time axis and satisfies: 

\\uj * (r) — u>o\\i < 5(e), as e — > 0, 

2 

where 5(e) is in Lemma 3.3. 

2. In the ball By(p) € H_ with p small enough and lo(0) £ By(p), there exists a stable 
manifold M_ in D(A?), codim.M_ = N, such that if the initial condition w(0) G M- 
and t — > oo, then 

\\u(t) - u*(t)\\ i < ce~ a l*l |K0) - w*(0)|| i . (3.15) 

2 2 

There also exists an unstable manifold A4 + in D(A^), dimA4 + = N, such that if the 
initial condition uj(0) € M.+ , then the inequality (3.15) holds as t — > — oo. In particular, 
if N = 0, then the unsteady quasi- geostrophic flow uj* is asymptotically stable. 

3. If the forcing f is almost periodic in D(A y ),^f > —\, then oj*{t) is almost periodic 
in D(Ai) with frequency basis contained in that of f . 

Proof. These assertions of the theorem follow from the representation 

uj = ujq + h(r, e) — ev(r, e) 

and Lemmas 3.1, 3.3. ■ 

This theorem gives stability estimate of stationary averaged QG flows; stability con- 
clusion of unsteady QG flows near the stationary averaged QG flow; comparison estimate 
between unsteady QG flows; and the existence of almost time periodic QG motions under 
almost time periodic wind forcing, on the entire time axis. 

Combining this theorem and earlier discussion in (1.15) and (1.16), we have 
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Corollary 3.5 Under the assumption of Theorem 3.4 and 



(3 2 \D\ 2 2\D\. 2 

4 ^>^^> n/oii<y^A?, 

2 \D\ 2 

where Ai = v — <-Ar~z as defined before, then for every R > such that the stationary 
averaged QG flow u)q has norm IIcjoIIi < R, the following inequality is valid: 

2 

\\u(t) -w*(t)||i < C(R)e- eaT \\u(0) -w*(0)||i,t -» oo, 

2 2 

where u>(t) = cj(t, e) is i/ie solution of the full quasi-geostrophic flow model (2.1) with 
initial condition uj(0) = ojq and e < €q(R). 

Combining Theorem 3.4 and Corollary 3.5, we conclude that there exists an a. p. 
solution for (2.1) when / is a. p., under the assumptions of Corollary 3.5. Here we give 
the restriction for ||/o||, not for ||/||. This result is obtained under different conditions 



from those of Duan and Kloeden [10] 



Let / be a.p. with values in H. Then by the results of |p5|- |^1 , also see (T(J, there 
exists a uniform attractor = ^h(/0)) °^ ^ ne non - au t onomous dynamical system (2.1), 
with TL = {f T ,f T (t) = f(t + r),r G \R}c b (\R;H)i an( ^ this attractor approaches to the 
attractor of the averaged dynamical system (2.3): 

Theorem 3.6 (Attractor convergence) Suppose that f is a.p. in H, then 

dist i (A H{fM) ,A) -> 0, asrj^oo, 

where A-j-i is the attractor for the quasi-geostrophic flow model (2.1), and A is the attractor 
of the averaged quasi-geostrophic flow model (2.3). 

This theorem claims that the uniform attractor of the full quasi-geostrophic flow model 
approaches to the attractor of the averaged quasi-geostrophic flow model, as the parameter 
in the oscillating forcing r] — > oo. 



4 Summary 

In this paper, we have discussed averaging principle for quasi-geostrophic motions under 
rapidly oscillating forcing (characterizied by a large dimensionless parameter r/), both on 
finite but large time intervals and on the entire time axis. We have derived comparison 
estimate, stability estimate and proved convergence result (as rj — > oo ) between quasi- 
geostrophic motions and its averaged motions. 

We also investigated the existence of almost periodic quasi-geostrophic motions under 
almost periodic forcing, and the convergence of the attractor of the quasi-geostrophic flow 
model to the attractor of the averaged quasi-geostrophic flow model as 77—5- oo. 
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